NASA (NGR) 33-010-189 No.S 


(NAS A-CR - 1 42 598) s ELASTIC HODDLI OF < 
PRECOHPRESSED PYROPHYLLITE USED IN ULTRAHIGH 
PRESSURE RESEARCH (Cornell lUniv.). ‘ 27 .p HC^ 
j 3 75 CSCL 20K 


G3/39 


N75- 20766 


Onclas 
18555. J 


ELASTIC MODULI OF PRECOMPRESSED PYROPHYLLITE 
USED IN ULTRAHIGH PRESSURE RESEARCH 


by 


Wolfgang Sachse 

Department of Theoretical and Applied Mechanics 


Arthur L. Ruoff 

Department of Materials Science and Engineering 

Cornell University 
Ithaca, N. Y. 14853 


December 1974 


Report #2366 
Issued by 



The Materials Science Center 



ELASTIC MODULI OF PRECOMPRESSED PYROPHYLLITE 


USED IN ULTRAHIGH PRESSURE RESEARCH 
Wolfgang Sachse 

Department of Theoretical and Applied Mechanics 
Arthur L. Ruoff 

Department of Materials Science and Engineering 

Cornell University 
Ithaca, N. Y. 14853 

ABSTRACT 

Chao and Ruoff [*] have shown previously that by precompressing pyro- 
phyllite prior to its use in high-pressure applications, one obtains a 
material which has a uniform density and uniform mechanical properties. We 
have recently studied the Dropaoation of ultrasonic poison in oyronhy 1 1 ir.p 
specimens to investigate the effect of specimen precompression on the measured 
elastic moduli. Measurements were made at room pressure and, for the precom- 
pressed specimens, to pressures of 3 kbar. 

We find pyrophyllite to be elastically anisotropic, apparently the result 
of the fabric present in our material. The room pressure adiabatic bulk modu- 
lus as measured on specimens made of isostatically compacted powdered pyrophyl- 
lite was determined to be 96.1 kbar. We found the wave speeds of ultrasonic 
pulses in pyrophyllite to decrease with increasing specimen precompression. A 
limiting value of precompression was found, above which no further decrease in 
wave speed is observed. For the shear wave speeds this occurs at 10 kbar while 
for the longitudinal wave at 25 kbar. In the limit, the shear waves propagate 
20 % slower than in the unprecompressed samples; for the longitudinal wave the 


[*] C. C. Chao and A. L. Ruoff, Rev. Sci . Instr., 4J_, 720-721 (1970). 



I. INTRODUCTION 


The use of pyrophyllite as an anvil - supporting material and as a 
pressure-transmitting medium in a Drickamer-type pressure apparatus is well 
known [1] Chao and Ruoff [2] have shown that by precompressing the pyrophyllite 
prior to its use, one obtains a material having uniform density and exhibiting 
uniform mechanical properties. They showed that this results in an improved 
work hardening response of the tungsten carbide pistons in the Drickamer pres- 
sure apparatus. In a recent publication [3], Ruoff has noted that the flow 
stress of many solids varies linearly with pressure, the rate involving the 
first pressure derivative of the effective elastic modulus. 

Although the characteristics and structure of pyrophyllite are well known 
[4-6], the authors have been unable to locate any publications concerning the 
elastic moduli of pyrophyllite or the pressure derivatives of these. Inasmuch 
as this precompressed material is expected to continue to find widespread 
application in ultrahigh pressure research, we studied the propagation of ultra' 
sonic pulses in pyrophyllite specimens, precompressed varying amounts to pres- 
sures of 40 kbar, to determine their elastic moduli at atmospheric pressure. 

Motivated by the aforementioned pressure-dependent flow-stress concept, 
we also measured the room-temperature pressure derivatives of the moduli in a 
range to 3 kbar above room pressure using 40 kbar precompressed pyrophyllite 
samples. 

II. EXPERIMENTAL 
a) Specimens 

Three types of pyrophyllite specimens having differing precompression 
histories were prepared. One series consisted of specimens as received, 
another consisted of specimens all precompressed to 40 kbar and a third con- 
sisted of a series of specimens each precompressed a varying amount which 



difference is 30%. 


The change in bulk modulus resulting from increasing pressure shows a 

behavior which can be approximated by a bi-1 inear curve with decreasing slope, 

the discontinuity occurring between 1.6 and 1.8 kbar. In the second portion 

of the curve, the bulk modulus derivative is determined to be 22.1. We believe 

that this value indicates the presence of voids which are closing under pres- 

1 /2 

sure. A regression analysis of the data shows that a P equation appears to 
give the best fit. The resulting empirical law for the bulk modulus derivative 
at pressure P becomes: 


dB = j L o 

dP [ (0.251P~0.147) 1/2 


(2.76P-2.56) 1/2 


J 


x 10 


-1 


where P is in kbars and in the present measurements , the longitudinal 

o 

modulus, is 204 kbar and , the average shear modulus, is 113 kbar. 

*o 
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ranged from 5 kbar to 40 kbar. PRECEDING PAGE BLANK NOT FILMED 

All specimens were cut from one portion of a large block of grade A pyro- 
phyllite purchased from the Maryland Lava Company. The specimens were fine- 
grained (smaller than 0.5mm) pyrophyllite having an observable fabric comprised 
of an alignment of the pyrophyllite grains and wispy greyish areas within the 
pyrophyllite. Visible under polarized light are several rectangular areas 
which may be pseudo-morphs after pyroxene because they exhibit mass extinction. 
Present also are small regions having limonitic stains, probably from the 
alteration of minute grains of pyrite, as well as spots of hematite coming from 
minute grains of magnetite. A micrograph of a typical section is shown in 
Figure 1 . 

The specimens comprising each of the three series listed above were cut 
sequentially from the length of the. original section. The dimensions of all the 
specimen discs were 2.233 ± 0.003 cm in diameter and 0.528 ± .003 cm in thick- 
ness. me dimensions of each specimen were measured with a Vani^euren iight-wave 

-5 

micrometer having a resolution of ± 1.3 x 10 cm. The faces of the specimens 

-5 

were parallel to ± 5 x 10 cm. 

b) Density Determination 

Most density determinations were made from specimen weight and dimensional 
measurements, but some were checked with a gravimetric method. 

Because pyrophyllite is so porous, we coated the specimen with a thin, 
water-proofing coating (Krylon), which, in all cases, represented less than 
0.05% of the specimen's weight when dry. We found that the amount of water 
absorption was limited to less than 0.17% of the specimen's weight after a five 
minute immersion in water containing a wetting agent. In all cases, the uncer- 
tainty in the density determination was less than 1%. 

c) Ultrasonic Measurements 

The velocities of pulses of ultrasonic longitudinal and shear waves were 
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measured using the pulse-echo technique. The radio-frequency pulse was gener- 
ated and received by a (Matec Model 6600) r.f. pulse generator and receiver. 
Coaxially plated, 0.127 cm-dia., 2MHZ X - cut and AC - cut quartz transducers 
were used to generate longitudinal and shear wave pulses, respectively. The 
bond was sodium salicylitate (sal ol ) . The display oscilloscope was used to 
make transit-time determinations with an accuracy of 1 to 2%. Higher accuracy 
measurements were made with the echo-overlap technique [7], Details of the 
system used for the present measurements have been described previously [8]. 

In this method, the overlap of any two echoes is achieved visually on the 
viewing oscilloscope by proper triggering of the oscilloscope and pulse gener- 
ator. The echo-overlapping procedure was repeated ten times from which an 
average transit time t and its standard deviation ± at were computed. Absolute 


accuracy of the transit-time determination in the measurements reported here 

is better than 1 part in 10^, while changes in pulse transit-time are accurate 

to 1 part in 10°. The rcoTii pressure adiabatic longitudinal and shear moduli 
s s 

denoted by and respectively were evaluated using Equation (1) 


C o 4 * C 0 * *0 


2 Z n \ 2 AP. 2a£„ 2at 

hr 1 ±( ^ + V 1 + t 2 * 

\ 0/ M 0 0 T o 


Here and ap q represent respectively the uncertainties in the initial 
wave propagation distance and p Q the initial specimen density. In an isotropic 
solid, C L corresponds to A + 2p and C s to p , where x is the Lame parameter and 
p is the shear modulus. 


The adiabatic bulk modulus B s and Young's modulus E s are given by the 


equations 

B S = <T - 4/3 C* (2) 

and 

E s = 3C|B S / (C* - Cj) 

Only changes in the pulse transit-time were measured as a function of 


( 3 ) 
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applied hydrostatic pressure. The corresponding change of elastic modulus is 
then determined when the pressure-dependent specimen dimensional changes are 
taken into consideration. Determination of this length correction can be made 
if the isothermal bulk modulus B is known. It is related directly to the adia- 
batic bulk modulus by the equation 

B S /B T = 1 + = 1 + A n ( 4 ) 

' pCp o 

The correction term a q contains 3 which represents the volume coefficient of 
thermal expansion (= 3a)(9.9 x 10” 6 ) [9]; B s is the adiabatic bulk modulus 

(96.1 kbar) [Appendix B] ; T is the absolute temperature (300°K); p , the 
density (2.738 g/cm^) and c p , which is the specific heat at constant pressure 
(0.22 cal/g°K) [10]. Using these, we find a q to be 1.12 x 10“^. 

The length at pressure P is then determined approximately for small length 
changes from 


P(l+A ) 

UtAf) = 1 + j — -f- dP (5) 

0 0 J o 3B 5 

so that the change in modulus including the length correction is determined 
from transit-time measurements according to 



Here t q and r represent the average round-trip time-intervals at atmospheric 
pressure and at pressure P , respectively. The uncertainty in this determina- 
tion of the change of the apparent modulus is given for our cylindrical speci- 
mens of initial dimensions r Q , l as 


t 2 
± 2 ( ~) 


Ar 


. A£_ At 

0_ + 0_ + 0 + At_ 

0 ^0 t ■ T 


( 7 ) 


We measured the attenuation of ultrasonic pulses travelling in pyro- 
phyllite to be quite high (15-18 db/cm). As the specimens were quite thin. 
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we found that 2MHZ pulses provided optimal material penetration and pulse 
definition so that accurate transit-time determinations could be made. 

The pressure system, using helium as a pressure medium has been described 
in detail previously [11], Pressure measurements were made with a manganin 
cell which had been calibrated against the new absolute pressure gauge system 
developed by Lincoln [12]. 

III. RESULTS 

A. Room Pressure Measurements (No Precompression) 

The average density of the pyrophyllite as determined from eight specimens 

o 

which had not undergone precompression was 2.675 ± 0.003 g/cm . From measure- 
ments of the longitudinal and shear wave speeds as a function of specimen loca- 
tion in the original pyrophyllite block, we found the average velocity of the 
longitudinal wave in all eight specimens to be 0.370 ± 0.004 cm/ysec. In each 
specimen we measured two distinct shear waves whose particle displacements were 
perpendicular to each other. The fast shear wave propagated with speed 
0.267 ± 0.003 cm/ysec while the slow shear wave propagated with speed 
0.233 ± 0.003 cm/psec. This indicates that the material is elastically aniso- 
tropic which is not surprising considering the pronounced fabric observed in 
the optical micrographs. 

B. Room Pressure Measurements (After Precompression) 

The average density of eight specimens which had been precompressed 40 kbar 

3 

was 2.738 ± 0.008 g/cm . From the wave speed measurements we found the average 
longitudinal wave speed to be 0.273 ± 0.004 cm/ysec and, as in the"as received" 
specimens, we observed two distinct shear waves which were measured to be 
0.215 ± 0.002 cm/ysec and 0.192 ± 0.002 cm/ysec. 

The pyrophyllite specimens which were precompressed varying amounts showed 
a permanent change in volume (increase in density) dependent on the 
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precompression load, similar to the results previously reported by Chao and 
Ruoff [2]. Shown in Figure 2 is the relationship between the volume ratio V/V Q 
and the precompression pressure - here V Q and V are respectively the volume of 
the disc prior to and subsequent to the precompression. Also shown in Figure 2 
is the change in ultrasonic wave speed in specimens precompressed differing 
amounts. The effect of precompression on a specimen is to decrease the measured 
longitudinal wave speed monotonically with increasing precompression. Pre- 
compression pressures in excess of 20-25 kbar were found to result in no addi- 
tional decrease in the wave speed. 

The intercept at zero precompression corresponds to that measured for the 
eight uncompressed samples, while the wave speed measured above 20-25 kbar 
agrees within 5 % of that measured for all eight 40-kbar precompressed specimens. 
The shear wave measurements show a similar behavior except that a preccmpression 
pressure of 10 kbar is sufficient to achieve the same shear wave speed behavior 
as observed In the 40 kbar precompressed specimens. 

These results, though contrary to empirical velocity - density relation- 
ships, c. f. Birch [13*14], are expected if structural changes result from the 
precompression of the pyrophyllite. A comparison of powder x-ray patterns of 
as received pyrophyllite and that which has been precompressed 40 kbar shows 
considerable line broadening of the latter, indicating that structural changes 
have, in fact, taken place. 

We investigated, in detail, the densification of pyrophyllite as a result 
of applied hydrostatic pressure in the range to 5 kbar. The specimens were in 
the form of cylinders 1 cm long and 1 cm in diameter. To insure against any 
exchange of gas between the specimen and the pressure medium, the specimens 
were compressed with a wet-bag technique in a pressure system in which mineral 
oil was the pressure medium. The specimen was held under pressure for thirty 
minutes to insure thermal equilibrium. After removal from the pressure vessel, 



-7- 


the specimen was unbagged, its axial deformation measured, rebagged and 
returned to the pressure vessel for the next loading cycle. Calculation of the 
volume change was made from measurement of the axial deformation and the 
assumption that our fine-grained specimens densified isotropically. The results 
obtained from three specimens is shown in Figure 3. Also shown in this figure 
are the volume ratios of the specimen precompressed 5 kbar and that precom- 
pressed 10 kbar, which were shown previously in Figure 2* 

We see that the change of volume appears to proceed in two stages, joining 
at approximately 1.8 kbar. These observations are similar to the predictions 
which Nelson [15] makes from consideration of a model describing the collapse 
of a single void in a plastically-deforming matrix. The first stage results 
from a non-linear stress-strain response of the elastically-deforming matrix 
while the second stage is predicted from consideration of the equilibrium 
between additional growth of the plastic zone of the matrix and collapse of the 

VOlu contained T n it. 

C. Pyrophyllite Modulus Measurements as a Function of Pressure 

Ultrasonic wave speed measurements could be made consistently to pressures 
of 3 kbar before the specimen-transducer bond failed. Measurements were made 
every 200 bars during pressurization and depressurization in 30 minute intervals 
to insure thermal equilibrium between the specimen and the pressure system. 

The changes of the apparent longitudinal and (fast) shear moduli, defined 
previously by equation (7), resulting from hydrostatic pressure are shown in 
Figure 4. Both curves exhibit an initial, linear change with pressure which 
decreases significantly at pressures exceeding 1 .6-1 .8 kbar. The drawn-in 
straight-line segments result from a linear regression analysis of the data in 
each region, above and below 1.8 kbar. 

The slopes of these lines are related directly to the first pressure 
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derivative of the elastic constant, being in fact, 

1 dC 
C o^ 

By assuming the shear modulus to be approximately equal to the average 
of the moduli specified by each of the two shear waves, we can evaluate the 
apparent Young's modulus E s and the apparent adiabatic bulk modulus B s for our 
specimens of pyrophyllite using equations (2) and (3). Our results are sum- 
marized in Table I. We see from Table I that the large change of modulus with 
pressure shown in Figure 4 results in an extraordinary large bulk modulus 
pressure derivative below 1.6 kbar. Above this pressure, that is, in the 
upper portions of the bi -linear curves, the derivative, though reduced, is 
still far greater than generally observed for most materials [16]. Bundy [17] 
determined B Q and from direct pressure vs. volume measurements of pyro- 
phyllite in the range to 20 kbar and he found B Q to be 267.2 kbar and B^ to be 
20.2 kbar. The B^ measured ul trasonical ly in the present investigation is far 
lower than that determined quasi-statically by Bundy. However, the pressure 
derivative B s measured in the range above 1.8 kbar is similar to that measured 
by Bundy. 

We believe that the low initial bulk modulus and its large pressure deriva- 
tive are the result of the presence of voids which are closing under pressure, 

Tt is possible that even the material which has been precompressed 40 kbar 
contains voids which open during the unloading of the precompression load. 

We attempted to coat a specimen with transducer attached using a silicone 
sealant (RTV) to prevent any possible exchange of gas between the specimen and 
the pressure medium. Unfortunately, the coating resulted in such a large in- 
crease in attenuation of the ultrasonic pulses which precluded all wave-speed 
measurements. 

As an alternative, we used a regression analysis with the measured modulus 
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data to determine an analytic form for the change of modulus of the specimens 
resulting from the applied pressure. A linear, parabolic, general power and 
exponential least-square regression analysis was made to determine which gives 
the most accurate description (as measured by the correlation coefficient) of 
the experimental data. Figure 5 shows the results of this analysis applied to 
the data presented in Figure 4. The parabolic equation in pressure P gives 
the best description which becomes: 


( C /C 0 -1) = -^ 1 + 


P/A 1 - 


/ 2 x 3" X 2 


1 


1/2 


where P is the pressure in kbars, and 


when C = C, 


when C = C, 


A-, = 6.26 


A 1 = 1.23 x 10* 


x 2 = 1.01 


a 2 = -1 .04 x 10 


_ n f on 
~ o . ULU 


Ag - i . i o 


Using the above and equation (2) we .can compute the adiabatic bulk modulus and 
its pressure derivative as a function of pressure P. We obtain for the latter 


C L 

di a _o S o x 10 -1 

dP (0.251 P-0.147 ) 1/2 (2. 76P-2. 56)V2 

This result is plotted to 50 kbar in Figure 6. It shows that the equation 
predicts the bulk modulus derivative to be less than 6.0 at pressures in excess 
of 30 kbar. We show for comparison the result which is obtained if the modu- 
lus change resulting from pressure is described by a general power law or a 
logarithmic law. 


CONCLUSIONS 


Our measurements clearly show pyrophyllite to be elastically anisotropic, 
apparently the result of the pronounced fabric present in our material. The 
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room pressure adiabatic bulk modulus as measured on an isostatical ly compacted 
specimen was determined to be 96.1 kbar. We found the wave speeds of ultra- 
sonic pulses in pyrophyllite to decrease with increasing specimen precompres- 
sion, A limiting value of precompression was found, above which no further 
decrease in wave speed is observed. For the shear wave speeds this occurs at 
10 kbar while for the longitudinal wave at 25 kbar. In the limit, the shear 
waves propagate 20% slower than in the uncompressed samples; for the longitu- 
dinal wave the difference is 30%. 

The change in modulus resulting from increasing pressure shows a behavior 
which can be approximated by a bi -linear curve with decreasing slope, the dis- 


continuity occuring between 1 ,6 and 1,8 kbar. In the second portion of the 
curve, the computed bulk modulus derivative is 22.1, We believe that this 
value indicates the presence of voids which are closing under pressure. A 
regression analysis of the data shows that a F equation appears to give the 


TIsr wnfiiUinn r.Mn-iy.'! a ru j a -f* ■» « ■Fpiv' 'f'hp .Ku-l U Tnodu 1 U c Howiwaf iuo OV'P- 
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diets values which are less than 6.0 at pressures in excess of 30 kbar. 


§ 



APPENDIX A: 


Orientation Dependence of Wave Speeds in Pyrophyllite 

The results of our shear wave measurements indicated that our material 
was elastically anisotropic. In an effort to determine what the elastic aniso- 
tropy is, we made additional wave speed measurements on specimens of unpre- 

i 

compressed pyrophyllite cut at various orientations to the discs described in 
Section II. To specify a rectangular coordinate system in the original pyro- 
phyllite block, we choose the axis parallel to the axis of the discs to be the 
x 3 axis. The x-j and x 2 axes are perpendicular to the x 3 axis and to each 
other. Plane-parallel specimens of approximately 0.5 cm thickness were cut 
normal to each of the three axes. 

The results of the wave speed measurements on each of the three specimens 
is listed in Table A. I. The results clearly show that the wave speed varies 
significantly with propagation direction. Furthermore, two distinct shear 


waves are observed in all three propagation direct ions. The shear waves sent 
in the x-j - direction exhibit their maximum responses in directions which are 
approximately 45° to the x^ and x^ directions. 

We attempted to determine from these measurements what matrix of elastic 
moduli would best describe the anistropic characteristics of the wave propaga- 
tion in our specimens. Because the longitudinal wave speeds in the x-j , x^ and 
x^ directions differ markedly, we conclude that the modulus matrix does not 
possess a cubic symmetry. 

By assuming the pyrophyllite to be transversely isotropic, we were able to 
determine the non-linear algebraic equations which specify the wave speed of 
the longitudinal or shear wave in terms of its propagation direction cosines 
and the elastic moduli. Substituting the measured wave speed data into these 
equations and knowing the geometric relationships between the various directions 
of wave propagation, we can determine the five elastic constants of our material 
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as well as the orientations of the specimens on which the wave speed measure- 
ments were made. We used a digital computer and Newton's Method to search for 
the roots of our equations. Unfortunately, the results never converged to 
meaningful values for the elastic constants or the direction cosines. This 
may be a consequence of the properties of the root-searching method employed 
or it may indicate that the data were not obtained from an elastically trans- 
versely isotropic material. Extensions of the procedure to higher elastic 
anisotropies were not made because it is clear that the computations become 
even more untractable. 

TABLE A. I 

Wave Speeds ± 0.010 (cm/ysec) 


Propagation 

Particle 

Displacement 

Direction 

Direction 

X 1 

x 2 

x 3 

x i 

0.306 

0.233 45° 

0.233 



0.215 to x^-x^ 

0.215 

*2 

0.231 

0.403 

0.256 

x 3 

0.218 

0.250 

0.345 



APPENDIX B 


Because we were unsuccessful in determining the elastic anisotropy of our 
material from the wave speed measurements, it was impossible to provide a 
method by which the adiabatic bulk modulus could be specified from these 
measurements. As an alternative, we attempted to produce an elastically 
isotropic sample of pyrophyl 1 ite. 

We produced a fine pyrophyl 1 ite powder (finer than 100 mesh) from the 
original bulk material by using a rock crusher. Compacting this in a wet-bag 
technique with an isostatic compaction pressure of 8 kbar, we obtained a 
"green" compact of pyrophyl! ite. Its strength was sufficient to permit the 
machining of several disc-like specimens of various orientations from the 
compact. 

Tht specimen's density was determined with the immersion technique which 
was described previously. The measured density was 2.576 ± 0.002 g/cm which 
represents 94.1% of the density of the pyrophyl 1 ite specimens which had been 
precompressed 40 kbars. From the wave speed measurements in various directions 
we found the longitudinal wave to travel at 0.361 ± 0.012 cm/ysec and the 
existence of one shear wave, travelling at 0.276 ± 0.010 cm/ysec, indicating 
that the material was elastically isotropic. 

We calculate from these data an effective longitudinal modulus (x + 2y) of 
335.7 kbar and an effective shear modulus (y) of 196.3 kbar. We applied the 
law of mixtures to the pyrophyl 1 ite-void specimen to determine the longitudinal 
modulus of a specimen if it consisted solely of pyrophyl 1 ite and its density 
was the same as that of the 40 kbar precompressed specimens. The adiabatic 
bulk modulus of such a material is calculated to be 96.1 kbar. 
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TABLE I 


0 < P < 1.6 kbar 1.8 < P < 3.0 kbar 


Modulus 

C(kbar) 

dC s 

jr 

Z-W< kh-' 1 ) 
0 

dC s 

W 

^(kbar' 1 ) 

c s 

U L 

204. ± 7.0 

40.39 

0.198 

36.20 

0.177 

c s 

113. ± 3.0 

8.85 

0.078 

6.97 

0.062 

B S 

53. ± 10.0 

28.59 

0.539 

22.13 

0.418 

E s 

198. ± 63.0 

53.54 

0.270 

31.44 

0.159 
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CAPTIONS 

1. Thin-section micrograph of Pyrophyllite. 

2. Relation between the room-pressure longitudinal and the shear wave 
speeds (2MHZ) as a function of pre-compression pressure of pyro- 
phyllite. Also shown (-Q-) is the change in volume ratio V/V Q as 
a function of the pre-compression pressure. 

3. Densifi cation of three specimens of pyrophyllite as measured by the 
room pressure volume ratio V/V Q resulting from a cyclic, increasing 
hydrostatic pressure. Also shown is the data of Figure 2 in the 
range to 10 kbar. 

4. Change of longitudinal and shear (fast) moduli as a function of 
pressure, measured during loading and unloading. 

5. Linear, exponential, parabolic and general power regression analysis 
of the data in Figure 4. 

6. Extrapolation of the analysis of Figure 5 to higher pressures. 
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Figure 1 
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Figure 4. 




Figure 5. 








